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Abstract. Wc have calculated the electrostatic interaction between two rod-like 
charged objects with arbitrary orientations in three dimensions, we obtained a closed 
form formula expressing the interaction energy in terms of the separation distance 
between the centers of the two rod-like objects, r, their lengths (denoted by 2li and 
2^2), and their relative orientations (indicated by 9 and (f>). When the objects have 
the same length {2li = 2I2 = I), for particular values of separations, i.e for r < 0.8/, 
two types of minimum are appeared in the interaction energy with respect to 9. By 
employing the closed form formula and introducing a scaled temperature t, we have 
also studied the thermodynamic properties of a one dimensional system of rod-like 
charged objects. For different separation distances, the dependence of the specific heat 
of the system to the scaled temperature has been studied. It is found that for r < 0.81, 
the specific heat has a maximum. 
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1. Introduction 

Electrostatic interactions between rod-like and stripe-like charged objects play a 
substantial role in the many systems of condensed matter and soft condensed matter 
physics such as strongly correlated materials, liquid crystals, electrolytes, polymers, and 
bio-molecules. For example, superconductors and anti-ferromagnetic Mott insulators 
are a vast class of strongly correlated systems which their quantum phases have been 
studied intensively. Several evidences signifying the segregation of charged carriers into 
fluctuating stripes in cuprate superconductors have been augmented [H [21 El S]. In 
the ceramic samples La2_a:Sr3.Cu04 co-doped with Nd the charged stripes are seen by 
neutron diffraction as a static lattice modulation coherent with the low-temperature 
tetragonal modulation [I1I2]- Nickelates P (La2_xSr2.Ni04) are the other good empirical 
realities for observing the charged stripes at low temperatures. The relationship 
between correlations of charged stripes and high-Tc superconductivity has been studied 
intensively in nickelates. More physical insights are gained by obtaining a closed form 
formula for the interaction of such stripes. Indeed, the closed form gives us many useful 
information on the physical properties of the system such as orderings and correlations 
of the charged stripes and also the thermodynamic properties of the system. 

Experimental evidences of such interactions are also observed in the colloidal world. 
Actually in the colloidal solutions, macromolecules such as DNA molecules, TMV or fd 
virus, V2O5 ribbons, crystalline Boehmite (AlOOH) rods, etc.[5l|6l[71[8] are intrinsically 
very anisotropic and have rod-like and ribbon-like shapes. In the context of orientational 
phase transition, it has been elucidated that the electrostatic interaction between the 
polyelectrolytes is lead to a twisting effect which enhances the concentration at the 
isotropic to nematic (I-N) phase transition P,[inj. More findings are attained by making 
use of the form of the interaction energy between such objects. Many efforts have been 
done to employ the analytical and numerical approaches for computing the electrostatic 
interaction of two rod-like charged objects. In a route developed work, D. Chapot 
et. al considered two anisotropic particles, separated by a large distance (i.e a distance 
larger than their typical dimension) and generalized the Derjaguin, Landau, Verwey and 
Overbeek (DLVO) calculations for the anisotropic molecules [12]. They have obtained 
an expression for the interaction energy between the particles [H]. However, because 
of performing various expansions, their expression is valid in the far field limit, i.e for 
inter-particles distance larger than both the Debye length (Ad) and the typical size of the 
colloid a. Meanwhile, the surface charge of the objects is involved in their calculations 
which should be obtained numerically. In the other hand, using molecular dynamics 
and Monte carlo simulations the interplay between the electrostatic interactions and 
orientational correlations have been studied for a ID array of charged rods [13]. 

Despite of all simulations and analytical approaches, the lake of a closed form 
formula for the interaction energy of such anisotropic particles is felt. In this paper we 
have obtained a closed form formula for the interaction energy of two rod-like charged 
objects with lengths 2li and 2I2 and arbitrary orientations in three dimension, separated 
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by r. By using the closed form formula we have computed the thermodynamic functions 
of a one dimensional array of charged rods in which each rod interacts only with two 
nearest neighbors. Although, the electrostatic interaction is long range and there are 
no practical realization to vindicate our assumption, one can get more insight about 
our finding closed form formula and its direct effects on the thermodynamic behaviors 
of a system. Finally the relationship between the behavior of the specific heat and the 
electrostatic interaction energy is investigated. 



2. Computational Details 

Let us consider a system of two inflexible rod-like charged objects with lengths 2/i and 
2/2, which are separated by a distance r. The interaction energy between two such 
uniformly charged objects is given by the following formula 



A1A2 e--^\^^-^^\dhdk 

47r£ J-iJ-i^ |r2-ri| 
where ri and r2 indicate the position vectors of two charge elements dqi = Xidli and 
dq2 = \2dl21 respectively. The inverse Debye screening length, k^, does characterize 
the decay rate of the interaction energy in the far field limit. This provides a way 
to find the screening factor from interaction force measurement. To calculate the 
interaction energy of two such rods it is convenient to work in the spherical coordinates 
(Fig.Q). Without lose of generality, we can set up the coordinate system in which 




Figure 1. Schematic representation of tlie orientations of two rods, separated by a 
distance r. The orientations of each rod are indicated by two types of angle. ^1(2) 
are the angles of the rods with z axis (co-latitude angles) and 01(2) are the angles of 
^1(2) sin 6*1(2) with x axis (azimuthal angles). 



the first rod is centered on the origin and aligned along the z axis and the second rod 
makes the angle 9 with z axis. Meanwhile the projection of the second rod on the xy 
plane makes the angle with x axis. 9 and are related to the old coordinates by 
cos 6* = cos 6*1 cos ^2 + sin^i sin6'2 cos(02 — 0i) and (j) = (p2 — 4>i- Indeed cos^ is ti • £2 
where ti and £2 are two unit vectors in the rods' directions and (ti • £2) = (cos 6*) shows 
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the correlation between the two rods' orientations. By introducing the polar coordinates 
for dli and dl2, M. Fixman and J. Skolnick calculated the screened coulomb interaction 
of two rigid rods of length /, where / assumed to be much larger than the Debye screening 
length [15]. Except for a proportionality constant, they put forward the result obtained 
previously by Onsager |TB]. In the limit of large Debye screening lengths, i.e — ^ 0, 
the screening effects are neglected and the exponential term in the interaction energy 
Eq.Q is close to one. In this case, by defining the dimensionless interaction energy, w, 
we can write the interaction energy as following: 



E, 



int. 



47re(2/i)(2/2 



-w. 



where. 



w 



E 

= l,2;a=± 



dzidz2 



-\Yz2VS^ ^\Z2 cos^ + a(-l)^2;i] 



(2) 



(3) 



Here we have written four kinds of integral in a compact form, by introducing the 
coefficients a's, the indices i's and dimensionless functions: H = sin ^ sin 0, J = 
sin 6^ cos 0, i\(2) = and dzi(^2) = ^?x^- Actually we have divided each rod to two 

1(2) 

equal halves. Then we have considered the interactions between each half a rod with the 
two halves of the other rod. Using some transformations, we calculated the integrals 
of Eq.([3]) analytically to attain to the closed form of the interaction energy. As an 
illustration let us describe the solution of the integral which is indicated by a = + and 
i = 1. By integrating out Zi the remaining expression is given by 

Z2 COS 9 + 



In 



l + zi + 2^2H 



In 



Z2 COS 9-£i + ^Jl + zl + il + 2z2H - 2£iZ2 cos 6 ) 



\dz2, 



(4) 



The first logarithm in Q is canceled by the same term with opposite sign in the other 
integrals which have not been presented here. However the second integral in Q is 
solved by part and is given by the following form 



In 



; COS 6* 



+ 



+ £2^ 



12 



dzo 



^+U^-ilf{z2) 



(5) 



f{z2) [z2COSe-£i + f{z2)] 

and X = 2(H-£icos0). Since (u^ - ^) > 



where f{z2) = + + Z2X, u 
mathematically we are allowed to perform the following transformation for calculating 
the integral of Eq. ([s]) 



Z2 + 



X 



(6) 



where y is an auxiliary variable. By employing the transformation ^ we can solve the 
integral of (|5|, simply (See the appendix). 

By making use of the same procedure to solve the other integrals in Eq. ^ i.e 
different a and z, we obtain a closed form formula for the interaction energy in terms of 
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the rods' lengths, their separations, and the angles 6 and 0, which specify the relative 
orientation of two rods. 

qiq2r 



int. 



47r£(2/i)(2/2 



w, 



1 



W = l2 

n=0 
1 2 



:Cos^ + £i + ^£j + £l + l + 24 [(-1)"H + ii COS 6 



COS e -£i + y/£l + il + l- 2^2 [(-1)"H + £i cos 6] 



n,m=0 j,k=l 



n+m 



H 



1 + (-l)'=cos^ 



In 



tanh 



sinh ^ 



+ (-l^^+^H 



-ly ^i cos^ 



-1)' 



+ 



2 cos 4 
sin^ 



E E(-i)" 



n,m=0 j=l 



1 

LJ 



P™ + (-1)^4 sin2^ + 



X 



tanh 



sinh ^ 



+ 



sm I 



sin^ 



E 



2 

Am,n tanh 



ni2 + (-1V'+'"H - {-lyii cose 



I sinh^ 



-lr m^2+(-l)"+'"H-(-l)"^icos6> i 

[ pm J 



+ P™COS^ 



PT^COS^) +J 



m,n=0 



Bm.„ tanh 



1 ginl^-l p^2-(-l)"+'"H-(-l)"£i cosg j 



P?,"COS0 +J2 



where 



A. 



m,n 



h sin^ e 



\n+m 



B. 



m,n 
m 



P™ + sin^^- (-1)"^ 



2' 
F 



(7) 



P™ = (1 -m)P+ + mP-, 
P™ = (1 -m)P" +mP+, 
P± = VG±4F, 
G = 1 -sin2^(sinV 
F = sin [2^] sin0. 

As it is explicitly seen from Eq. ([T]) the closed form formula depends on the variables r, 
9 and (f). This formula is utilized as the starting point for the most computations in the 
various subjects of condensed matter and soft condensed matter physics such as the field 
of strongly correlated materials, liquid crystals, polymers and biophysics. Employing 
the closed form formula ([T]) we can investigate the thermodynamic properties, the order 
parameters and the different phases of a system, more precisely. It should be noticed 
that the above expression for the interaction energy of two rod-like particles is realistic 
for separations larger than the particles size, i.e the separations larger than the diameter 
or lateral extension of the particles. 
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Let us scrutinize the different rods' orientations with respect to their interaction energy. 
In Fig. ([2]) the scaled interaction energy have been plotted versus 9 and different values 
of (p and r. Let us focus on the solid line curve corresponding to = | and r = 0.21 
in Fig.([2]-a). Considering n as an integer number, it is clearly seen that the interaction 
energy for the case of 6 = rnr is stronger than for the case of 6 ^ nn and for the case of 
9 = ^^^Y^'Tc is weaker than the case of 7^ ^^l^^- Moreover, the configurations indicated 
by {9 = mr, (j) = = 0.21) and {9 = ^^^vr, = |, r = 0.21) are the maximum and the 
minimum energy configurations, respectively. By increasing the interaction energy 
is heightened, except for the directions of 9 = rni. Indeed the situation denoted by 
9 = nil corresponds to the parallel rods and does not depend on 0. Thus no variation 
is seen on the energy value of the two rods and w{9 = mr, 4>,r = 0.21) ~ 15. For the 
orientations which are indicated by ^ 7^ nir, an enhancement of makes one side of the 
tilted rod closer to the fixed rod and the interaction energy is increased. The appearing 
of two minimums, viz. the local and the global minimums in the interaction energy 
are noticeable (see FigMa-dashed line). Therefore by increasing two minimums are 




e e 

Figure 2. The interaction energy between two rods with equal lengths 2li = 2I2 = 
I — 8 versus 9 and different values of (j) a-nd r. 

emerged on the interaction energy at = vr (when two rods are parallel) and ^ = | 
(when tow rods have been lied on two individual perpendicular planes). To see the 
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configurations with the local or the global minimum energy, we have divided the whole 
region of </> G [0, vr] to the five intervals (see the table ([T])). 



Table 1. The minimum and the maximum energy orientations for the rods which have 
been separated by r = 0.2L G-Min: global minimum and L-Min: local minimum. The 
G-Min value is smaller than the L-Min one. 








/ 7r TZ \ 


1-2 ' 8 / 


r57r Stti 
I 8 ' 4 J 


( — >7^J 


6 = m: 


Max. 


L-Min. 


G-Min. 


L-Min. 


Max. 


0= 2"+i7r 


Min. 


G-Min. 


L-Min. 


G-Min. 


Min. 



Let us compare the interaction energy between two parallel and two perpendicular 
rods which are separated by the distance r = 0.21. As it is obviously seen from the table 
([l]) at the interval (j) G [|, ^) the parallel orientation is the global minimum energy 
configuration, i.e the rods prefer to be parallel, however in the other intervals of </> the 
global minimum energy is gained when two such rods are perpendicular. These intervals 
are not universal and depend on the distance of the separations. To get more insightful 
information we have also plotted in Fig.(|2]-6, c and d) the scaled energy w versus 6 for 
different separation distances r. For r > 0.7/, the orientations indicated by {6 = rni^cf)) 
and {9 = '^^y^tt, (p) are always the minimum and the maximum energy configurations, 
respectively. Thus for r > 0.7/ two rods prefer to be parallel. By decreasing the distance 
of the separation from 0.7/ the relative orientations corresponding to the maximum and 
the minimum energies are changed. 




Figure 3. $im: The interval of (j) where the orientations {9 = mr, $im) are the global 
minimum and ^2m- The interval where two minimums are appeared in the interacting 
energy. 

To discover the minimum energy orientations it is more beneficial to introduce an 
interval, say $im in which the orientations indicated hj 6 = nir (parallel orientations) 
are the global minimums. We have plotted the function $1^,/ versus r in Fig.(|3]). As 
it is clearly seen is increased by increasing r and saturated approximately at 
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r ~ 0.8/. It is also useful to define another interval, $2m, where the state of the system 
is characterized by two minimums in the interaction energy. The corresponding relative 
orientations are indicated by {9 = nir, $2m) and {9 = ^^y^vr, $2m)- As it is also seen from 
the Fig. ([3]) by decreasing the separated distance r, the intervals of $2m are increased 
and reached to a maximum value at r ~ 0.3/. It is remarkable that the interval <l>2m has 
considerable values only for r < 0.8/. This characteristic behavior gives an insightful 
intuition when we discuss the response function of a ID array of the interacting rods. 

4 
3.5 

3 
2.5 

1.5 
1 

0.5 

4 8 12 16 20 

r 

Figure 4. The interaction energy between the rods {I = 8) versus r and different 
relative orientations. 

Fig.(|4]) shows the interaction energy of two charged rods as a function of their 
separations r, and different relative orientations. The interaction energy decreases with 
increasing the distance of separation between the rods. It is seen that all curves coincide 
at large r values, indicating that at large values of r, the rods can be treated as point 
charges, which is in well agreement with other works [llj. 

To see the advantage of the closed form formula, let us consider a ID array of 
interacting rods which are separated by a lattice constant r. Typically this model 
is used to discuss the properties of the ID liquid crystals. We have considered the 
situation in which each rod interacts only with its nearest neighbors (N.N). Indeed, the 
electrostatic interaction is long-range and there are no practical realization to vindicate 
our assumption. However we bring this example, to show the direct influence of the 
closed form formula ([T]) on the thermodynamic behaviors of the system. 

We have defined a re-scaled temperature, t = qf^ , where T is the temperature, 

fee is the Boltzmann constant, and £ = Ms a dimensionless length (/ = 2/i = 2/2). The 
re-scaled temperature can be expressed in terms of Bjerrum length is = e"^ /AnekBT 
(For water £b = 7A at room temperature). It can be seen that t = ^ (f ) , where n 
is the number of fundamental charges (e) on the rods. Fig.(|5]) shows the change in the 
specific heat per particle as a function of the re-scaled temperature, for a ID array of 
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Figure 5. The specific heat for a chain of N.N interacting rods (/ = 8). Right: for 
r < 0.7/ where two types of minimum are appeared in the interacting energy plot. 
Left: for r > 0.71 



N.N. interacting rods. 

Let us construe the behavior of the specific heat with respect to the interaction 
energy. As it is obviously seen from Fig. (|5|, For distances r < 0.71 there is a peak 
on the specific heat at low temperatures. However, this is not the case of r > 0.7/ and 
the peak disappears at low temperatures. Indeed this peak is seen when there exist two 
types of minimum in the interaction energy. Moreover, the width of the peak depends 
inversely to the value of ^2m, i-e at large values of ^2m, such a peak is very broaden and 
becomes narrower by decreasing the interval $2m- 

For more sense, it is worth to discuss the peak values and peak locations in 
the actual units. As it is observed from the specific heat data, for the separation 
r = 0.4/ {£ = 2.5) the peak value and peak location are c ~ 1.02 ks and t ~ 0.06, 
respectively. For a water solution, and highly charged particles ^ ~ 0.5 and the peak 
location corresponds to the temperature T ^ 10 K. It means that at T ^ 10 K the 
specific heat for a ID array of N.N rods has its maximum value. In the case of ^ ~ 1 
the peak occurs at T = 40 K. 

It is remarkable again that for a realistic model we should consider the interaction 
between more neighbors. Discussing on the thermodynamic behaviors of such models 
by using our closed form formula is our future work. 



Appendix 



Using the transformation ([g]), we can find the following expression for the integral of 
Eq.g: 



w = £2 1 - In 



2 cos^ - £1 + Jil + £i + 1 - 2^2 [H + h cosO] 
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1 

n=0 

2 cose 
smfc' 



H 



1 + (-1)" cos 61 



(-1)% 



In 



tanh [| sinh ^ [ 



fe+H-^i cos< 
P+ 



■]]-(-ir 



1 

J](-l)"tan-i -(Atanh 

n=0 



tanh [|sinh-^ [ H-^i^cosg jj _ (^_^^^ 
"n4 + H-^i cos 6' 



sinh ^ 



P+ 



+ p+ cos e 



)] 



sm , 
+ ^ln 



(A tanh sinh-i [ "-^pV'°'1 ] + P+ cos ^) V 
(A tanh [I sinh"^ [ ^2+H-|icose ]j ^ p+ ^) ^ js 



(8) 



where ^ = P+ + £i sin^ 61 + 1. 
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